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We propose a hybrid (continuous-discrete variable) quantum repeater protocol for distribution of 
entanglement over long distances. Starting from entangled states created by means of single-photon 
detection, we show how entangled coherent state superpositions, also known as 'Schrodinger cat 
states', can be generated by means of homodyne detection of light. We show that near-deterministic 
entanglement swapping with such states is possible using only linear optics and homodyne detectors, 
and we evaluate the performance of our protocol combining these elements. 



Two regimes for optical quantum conTputation and 
communication are commonly identified P, In the 
discrete variable regime, single photons are the carriers 
of information with qubits encoded e.g. in their polari- 
sation or frequency, and measurements are performed by 
single-photon detection (SPD). In the continuous vari- 
able regime, information is encoded in continuous degrees 
of freedom of the electromagnetic field such as the field 
quadratures, which are measured via homodyne detec- 
tion. Both regimes have their advantages and drawbacks. 
On the one hand, single photons are good for heralding 
events across lossy channels since their discrete nature 
implies that no partial loss can take place. Successful 
detection of a photon at the channel output unambigu- 
ously identifies a successful event (in the absence of dark 
counts). This is useful for entanglement generation 
On the other hand, certain tasks, such as quantum tele- 
portation [5], can be accomplished unconditionally with 
linear optics in tlie continuous regime while not in the 
discrete one, and homodyne detection allows quadrature 
measurements with much higher efficiency than what can 
be achieved with SPD at present. 

In quantum repeater schemes, entanglement distribu- 
tion across a long lossy channel is achieved by combining 
entanglement generation over shorter segments with en- 
tanglement swapping and purification [6|. Reaching high 
communication rates thus requires both a good genera- 
tion scheme and efficient swapping. Repeaters combining 
continuous and discrete variables have recently been pro- 
posed. Ref. 01 relies on non-linear interactions of light 
with single spins in cavities, which is not easy to realise 
experimentally. Ref. Q is based on linear optics, but 
makes use of inefficient single-photon counters to perform 
entanglement swapping. Here we describe a quantum re- 
peater (Fig. [1]) which combines entanglement generation 
based on SPD with continuous variable techniques to ef- 
ficiently distribute entanglement in the form of states re- 
sembling 

Wo{a))=e^'\a)\a)+e-''\-a)\-a}, (1) 

where \a) is a coherent state and is a phase (through- 
out this paper a is real). In quantimr optics, coherent 




FIG. 1: I Generation of Bell-type entanglement between two 
nodes. Light from two-mode squeezing sources is mixed on a 
balanced beam splitter and a single photon is detected. The 
remaining modes are stored in quantum memories (QM). II 
Generation of approximate two- mode cat states. Pairs of Bell- 
like states are joined at balanced beam splitters, and the X- 
quadrature is measured at one output of each beam splitter. 
When the sum of the outcomes is close to zero, the state is 
kept, and the process is iterated. Ill (a) Simple entangle- 
ment swapping. Two two-mode cat states are joined on a bal- 
anced beam splitter, and the outputs are homodyned. Success 
is conditioned on an X-outcome close to zero. Ill (b) Im- 
proved entanglement swapping using k auxilliary single-mode 
cat states of size (left to right) 2a,. . . ,2^1'^ a. 

state superpositions of this form, and similarly for a sin- 
gle mode 

|C'(a)) = |a) + |-a), (2) 

are often called 'Schrodinger cat' states. Such states are 
useful for a number of applications in quantum informa- 
tion, including fault-tolerant linear optical computation 
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and quantum teleportation ■ Generating exact cat 

states is difficult, since it requires e.g. very strong Kerr 
non-linearities Approximate cat states can be gen- 
erated by photon subtraction, as demonstrated recently 
by several groups [l^ - [l5| . but unfortunately the average 
photon number is usually restricted to < 1. 

We demonstrate that near-deterministic entanglement 
swapping of cat states using only linear optics and homo- 
dyning is possible, and we devise an efficient probabilis- 
tic scheme for generating states with many photons and 
very good overlap with exact squeezed single- and two- 
mode cat states. Our generation scheme is reminiscent of 
that of Ref. fli'l , where creation of squeezed single-mode 
cat states from photon-number states was demonstrated, 
but does not require input states with more than a sin- 
gle photon in each mode and takes advantage of quan- 
tum memories to significantly increase the rate. The ele- 
ments of our repeater protocol are outlined in Fig. [TJ In 
step (I), Fock state entanglement between two separated 
nodes is generated by means of SPD. In step (II), beam 
splitters and homodyne measurements are used to cre- 
ate two-mode cat states from the states generated in (I) , 
and in step (III) , the entangled cat states are swapped to 
larger distances. We discuss these steps one by one and 
then finally consider the protocol in its entirety. 

As shown in Fig. [TJ the entanglement generation step 
(I) can be implemented using two sources of two-mode 
squeezed states, realised using e.g. parametric down con- 
version crystals or ensembles of A-type atoms 0, [3| ■ 
One photonic mode from each source impinges on a cen- 
tral beam splitter, and a single photon is collected as 
indicated. Following a single SPD click, the remaining 
two modes are projected to a Bell-like state 01 
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(|oi) + |io)) + o(v?), 



(3) 



where the last term represents contributions from mul- 
tiple excitations and is small when the pair production 
probability p is small, i.e. for weak squeezing. A setup 
of this type benefits from the discrete nature of single 
photons, with clicks heralding successful transmission. 

In step (II), states of the form ([3]) are combined on 
balanced beam splitters, and the A"-quadratures of two 
of the output modes are measured. The resulting state 
is kept, if the outcomes sum to a value within an interval 
[—A, A]. Otherwise it is discarded, and the process is 
restarted. Upon success, the protocol is iterated with 
the output states as new input states. To see that this 
scheme can produce cat-like states, it is illuminating to 
first consider the limit A — 0, p — >■ and generation 
of a single-mode state (illustrated by, say, the left hand 
side of Fig. [l]II only). We start from two sources, each 
producing a single excitation |1) corresponding to the 
standard harmonic oscillator wavefunction 



The joint wavefunction for both sources has the form 
ipo{x)'4'o{y)- A balanced beam splitter (any non-diagonal 
unitary will do) is then applied to the pair of modes 
x,y, transforming the state to ^poii^ + y) / \/2)'4'o{{x — 
y)/y/2) oc e~^^ +^ ^/^(a;^ — y^), followed by a measure- 
ment of y. If we require y — 0, corresponding to A — >■ 0, 
the resulting unnormalised output state takes the form 
e~2^ x'^. The process is now iterated, combining this 
state with the output from another pair of sources, etc. 
After m iterations (in Fig. [T]II, m = 2), the final nor- 
malised output state wavefunction becomes 



V'™(a;) =r(2" + l/2) 



^1/2, 



(5) 



(4) 



This expression is a symmetric, double-peaked function 
of X, which closely resembles the wavefunction of the 
state in Eqn. In fact, defining |^(a)) to be the (nor- 
malised) single-mode cat state, we find that |V'm) is well 
approximated by S'(2)|^(/i„i)), where /!,„ = ^2™ -1-1/2 
and S{s) denotes squeezing in the variance of X by a fac- 
tor of s. The output state is thus very nearly a squeezed 
cat state. The fidelity |(V'm|'S'(2)|^(/i,„))p of the actual 
output with respect to the squeezed cat exceeds 99% for 
TO > 2. Being squeezed does not render our states less 
useful. They are squeezed in amplitude, making them 
more robust against decoherence [19], and if desired they 
can be unsqueezed by local operations. The required 
squeezing of about 3dB or less is easily accessible in ex- 
periment. If the squeezing is introduced into the input 
Fock states, the single-mode setup corresponds to am- 
plification of small cat states via homodyne detection 
|20l |21 |. The scheme can be generalised from a single 
to an arbitrary number of modes by taking the mode of 
the single excitations from the sources to be a superpo- 
sition of some other set of modes. Eq. ^ for example 
corresponds to an excitation in a symmetric superposi- 
tion of two spatial modes. More generally the new mode 
variables form a vector x, and the source mode is given 
by a^x where a is a unit vector. For a pair of sources, 
the single-mode scheme is applied to each pair of cor- 
responding modes and we require |a''^y| < A, where y 
is the vector of measurement outcomes. Separating out 
the mode determined by a, the final state takes the form 
\i'm)a\vac), where \vac) is the vacuum state of the re- 
maining modes. 

Step (II) may be implemented entirely on traveling 
light beams (see Fig. [1]), in which case simultaneous suc- 
cess for all measurements is required. The rate can, 
however, be significantly increased if the quantum states 
generated at each level of the protocol can be stored in 
memories, since in this case there is no requirement for 
simultaneity. The required homodyne measurements can 
be performed on retrieved light fields or via quantum 
nondemolition measurements directly on the memories 
(2^ . To get a non-zero probability for successful genera- 
tion one needs to take a non-zero value of A, and there 
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FIG. 2: (Colour online) Probabilistic generation of approxi- 
mate cat states. The fidelity-rate trade-off is plotted for per- 
fect input states (dots) and input states with a 1% two-photon 
contribution (circles). The rate is measured in units of the 
source repetition rate. 

will be a trade-ofF between the rate of generation and 
the output state fidelity. To illustrate the trade-off with 
memories, we show in Fig. [5] a numerical simulation of a 
single-mode-cat generation. In the simulation, we include 
a double excitation contribution of 1%, and we find that 
states with a fidelity of 90% with respect to S'(2)|^(/i3)), 
corresponding to an unsqueezed cat amplitude a = 2.9, 
can be generated at a rate of ~ 0.08 measured in units of 
the source repetition rate, i.e. the rate at which (H)) can 
be produced. 

As illustrated in Fig. [T] Ilia, once entanglement has 
been established within two neighbouring segments, they 
can be connected by locally mixing the modes at the 
central node on a balanced beam splitter and measuring 
the X- and P-quadratures at the two output ports. We 
will first consider entanglement swapping with ideal cat 
states and subsequently with the actual states produced 
in step (II) of our scheme. We denote the normalised 
equivalent of |7g(a)) in Eq. ([1]) by \^e{a)). Note that 
the amount of entanglement in |7g(a)) varies with 9, but 
for coherent state amplitudes |ap > 2 the variation can 
be safely neglected and |7e(a)) contains one ebit of en- 
tanglement. We consider entanglement swapping of two 
copies of |7o(q!)). The idea of the swapping procedure 
is that several terms in the input state may lead to the 
same measurement outcome. Conditioning on this out- 
come thus projects the output modes onto a superposi- 
tion state, which is entangled. For coherent state inputs 
|ai),|a2)j the balanced beam splitter outputs coherent 
states with amplitudes (ai ±a2)/\/2- Since each input is 
I ± a), the possible outputs are | ± v^a) and |0). There 
are two ways to obtain the latter, and hence measuring 
X and conditioning on an outcome close to zero results in 
an entangled output state. More formally, the state prior 
to swapping is |7o(a))|7o(a)). Introducing the shorthand 
|a+) = |Q;)|a), |a+) = |Qf)| — a) etc., the (unnormalised) 
state after a P-measurement with an outcome po becomes 

\a^_)\V2a), + \a-)\ - V2a). + |7eo («)) |0)., (6) 



where we have labelled the mode for which X is mea- 
sured by X, and Oq = — 2apo- If a is large enough for 
|0) and \V^a) to be nearly orthogonal, we see from © 
that an X-measurement with an outcome close to zero 
will project the remaining modes to a perfect two-mode 
cat state. Since {7'gg{ct)\7gg{a)) ~ 2, the probability for 
this successful outcome is roughly 1/2. 

Probabilistic entanglement swapping of two-mode cat 
states is thus easy to implement. As we now demon- 
strate, the swapping can be made nearly deterministic 
using auxiliary single-mode cat states as a resource. The 
setup is shown in Fig. [T]IIIb. Additional beam splitters 
are inserted between the first beam splitter output and 
the X-measurement. At the j'th beam splitter an auxil- 
iary cat state |^(2^/^a)) is injected, and a P-measurement 
with outcome pj is performed in one output port. We can 
understand what happens by considering just one auxil- 
iary state. By mixing the output from the first beam 
splitter with a single-mode cat of amplitude \/2a we ar- 
range that each of the 'failure' outputs above can combine 
with one term of the cat state to yield |0). An X-outcome 
of zero then projects the output to an entangled state. 
At the same time, the 'success' output above splits into 
I ± a), leading to two possible measurement outcomes 
each of which still produces the entangled state loggia)). 
We have thus increased the number of X-outcomes which 
lead to successful swapping. Formally, after the two P- 
measurements the (unnormalised) state is 

l7^,(«))|0>. + he„{a)){e^''\a)^+e-n " (7) 
-t-|a+)|2a), + |a;)|-2a)„ 

where 9i — — 2'^/^api, v is an unimportant phase and 
\%^{a)) — e^^\at)+e^''^\a^) equals 17^^ (a)) up to a local 
phase shift. Assuming that |0) and \a) are nearly orthog- 
onal, X-outcomes originating in the |0) or | ± a) terms of 
([7|) all project the output to a two-mode cat state. Only 
the extremal outcomes lead to failed swapping. Counting 
terms, the success probability is seen to be 3/4. Gener- 
alising this to k auxiliary states, we find that the success 
probability scales as 1 — 2~''~^. Fixed distinguishabil- 
ity of the terms requires a to scale as a ^ 2^^/^. The 
failure probability thus scales inversely with the mean 
photon number in the two-mode cats and the square 
root of the mean photon number in the largest single- 
mode cat. This result demonstrates that, given suffi- 
ciently large cat-state resources, entanglement swapping 
using only linear optics and homodyne measurements can 
be performed with success probability arbitrarily close to 
one, i.e. near-deterministically. 

We now turn to the assembly of all three steps - gener- 
ation of entanglement in steps (I) , (II) and entanglement 
swapping with the generated states in step (III) ~ into 
a complete quantum repeater protocol. The goal of the 
protocol is to establish a useful entangled state across 
some channel of length L. To this end, the channel is 
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FIG. 3: Rate (solid line) of the full repeater protocol, includ- 
ing both generation and connection of entangled states. We 
assume a fixed final fidelity of 90%, a fiber attenuation lengtii 
of 20km and an SPD efficiency of 50%. The rhs axis shows 
the optimal values of n and m. 



divided into 2" segments of a shorter length Lq- En- 
tanglement is generated in each segment separately and 
they are then connected by entanglement swapping. We 
will consider only simple swapping (Fig. [T] Ilia) with no 
auxiliary states. The benchmark for the performance of 
a repeater is the rate at which final, entangled states 
can be generated over a distance L with a fixed fidelity 
with respect to some ideal target state. Just as in step 
(II), one must choose a finite acceptance range [—(5,(5] 
for the X-measurements during swapping, and for fixed 
L there will be a trade-off between repeater rate and fi- 
delity through the parameters A, 5, p. The target state 
is approached as A,(5, p — > 0. In this limit, the entan- 
gled states before swapping are given by \'iprn) +\vac) ^ . 
where ± stand for the symmetric and antisymmetric com- 
binations a ± & of the spatial modes a, b. Conditioning 
on zero in all X-measurements, we find that for suffi- 
ciently large m, up to two P-space displacements which 
can be trivially canceled, the wavefunction after n con- 
nections takes the form S+{^)S-{^)\^ip^ (Mm/ V^kn))- 

Here k„ = 2A/2coth(2"arccoth(l/-\/2)) while ipn depends 
on the P-measurement outcomes. Because fc„ converges 
fast towards 2^/2, this state is essentially equivalent to 
SaiV2)SbiV2)\j^,X^im/2^^'^)) . Our target state is thus a 
squeezed, nonlocal cat state. We remark that the 1.5 dB 
squeezing in the final state is local and can be locally 
undone, even though the squeezed mode for the states 
generated in step II is the nonlocal, symmetric mode '-I-'. 

We numerically simulate the repeater proto- 
col, optimising the rate £IS £L function of distance 
for a fixed final fidelity of 90% with respect to 
SaiV2)SbiV2)\-f-^„{tim/2"''/'^)). We optimise over the 
acceptance intervals A and (5, the number of generation 
and connection steps m, n, and the pair production 
probability p. For runtime reasons we restrict the 
simulation to m < 3. We take into account losses during 
entanglement generation, assuming a fiber attenuation 
length of 20km and an SPD inefficiency of 50%, and 



we take the source repetition time to be given by the 
classical communication time Lq / c, where c is the speed 
of light. The result is shown in Fig. O At 1000km 
our protocol reaches a rate of 0.3 pairs/minute. The 
best proposed atomic-ensemble-based repeaters in the 
discrete variable regime achieve comparable rates. 
However to do so they require either high efficiency 
SPDs (> 90%), which unlike the efficient homodyne 
detectors employed here are not readily available in the 
lab, or more involved entanglement swapping procedures 
[2^ . 24 1 . In addition, for repeaters based on linear optics 



and SPD in the discrete variable regime, the success 
probability for entanglement swapping can never exceed 
1/2, whereas for the protocol presented here we have 
demonstrated that near-deterministic swapping and 
hence a much higher rate is possible when the swapping 
procedure of Fig. [T] Ilia is replaced by that of Illb. An 
analogous setup can also be used for near-deterministic 
teleportation. Thus the hybrid approach represents a 
promising avenue for reaching higher rates for entangle- 
ment distribution. Since our protocol uses only linear 
optics, light storage and retrieval, and single-photon and 
homodyne detectors, the means for a first experimental 
implementation are available in today's laboratories. An 
interesting variation of the present scheme might be ob- 
tained by integrating it with the proposal for generation 
of two- mode from single-mode cats of Ref. Q . 
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